Possible Lattice Distortions in the Hubbard Model for Graphen^ 
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The Hubbard model on the honeycomb lattice is a well known model for graphene. Equally 
well known is the Peierls type of instability of the lattice bond lengths. In the context of these 
two approximations we ask and answer the question of the possible lattice distortions for graphene 
in zero magnetic field. The answer is that in the thermodynamic limit only periodic, refiection- 
symmetric distortions are allowed and these have at most 6 atoms per unit cell as compared to two 
atoms for the undistorted lattice. 

PACS numbers: 71.10.Fd, 61.48.Gh, 73.22.Pr, 05.30.Fk 



Graphene is a two-dimensional array of carbon atoms 
arranged in a honeycomb lattice. The number of delo- 
calized electrons, N, equals the number of carbon atoms, 
i.e., the conduction band is half- filled. The dynamics of 
the electrons is often modelled by the Hubbard Hamilto- 
nian 

HiT) =- J2 T.^-V i^l'rCy.a + h.C.) (1) 



The Cx^a are fermion annihilation operators and n^.o- = 
c\. a-Cx,a- For the purposes of this paper the on-site repul- 
sion U can have any sign as long as it is the same for all 
lattice sites. In particular, for = we have the Hiickel 
model. 

The notation emphasizes the dependence on the near- 
est-neighbor hopping matrix T = {txy)- We take the 
hopping matrix elements txy to be positive, although not 
necessarily independent of the pair x, y. (By the well- 
known particle-hole transformation on a bipartite lattice, 
we could, as well, take the t^y's negative.) It is important 
for us that txy is real, i.e., that there is no magnetic field. 

The interesting quantity, which does not usually ap- 
pear in the Hubbard model, is the distortion energy 
Fitxy). What we are assuming is that txy depends on 
the physical distance dxy between the lattice sites x and 
y. There is an equilibrium value and deviations from 
this value cost a positive energy. After eliminating dxy 
we can assume that the distortion energy depends on 
txy This is the quantity F{txy) in Thus, F{t) has a 
minimum ai t — t'^^\ corresponding to d^^\ and it goes 
to infinity as t goes to zero (infinitely separated atoms) 
or to infinity (no atomic separation). We do not have 
to assume that F is convex. This model, while quite 
general, does not take account of the possible distortion 
energy connected with the change in bond angles or pos- 
sible non-planar atomic configurations, i.e., curvature of 
the lattice. It also does not take account of longer range 
interactions 0, [I]- 

The problem we address is the possible configurations 
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FIG. 1: (a) On the left is the only graphene structure allowed 
by our theorem. The locations of the three possible bond 
lengths are indicated by heavy-solid, light-solid and dashed 
lines. There are six atoms (or three hexagons) per unit cell. 
The Kekule structure is the special case of equal length solid 
lines, (b) The ALT structure on the right breaks some of the 
reflection symmetry and is not allowed by our theorem. 



T = (txy) that minimize the ground state energy E{T) 
of H{T). We answer this in the thermodynamic limit by 
proving that one ends up with a relatively simple peri- 
odic configuration. Nothing more complicated can occur! 
This periodic configuration is shown in Fig. 1(a), where 
it will be seen that the resulting lattice can have at most 
three different lengths and the unit cell can have at most 
six atoms. A special case is the so-called Kekule lattice 
in which two of the lengths are equal. It also has six 
atoms per unit cell. Many authors have considered vari- 
ous possible energy-minimizing periodic distortions. The 
discussion can and should also be extended to include 
non-periodic and chaotic structures. Our contribution 
is to rigorously exclude anything more complicated than 
that shown in Fig. 1(a), in the thermodynamic limit. In 
particular, the ALT structure of Ref. [3|-|5| in Fig. 1(b) 
is not an energy-minimizer for the Hubbard model, since 
it breaks two of the three reflection symmetries shown in 
Fig. 1(a). 

The idea that distortion of a lattice can lower the en- 
ergy is sometimes called the Peierls instability, but the 
idea is older (see @). Peierls' essential contribution 
was to realize that this distortion always occurs in ID 
for an infinitely long chain or ring of atoms. It does not 
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always occur in higher dimensions, or in a finite ID ring, 
e.g., benzene. 

Let us recaU the situation in one dimension, i.e., the an- 
nulenes hke benzene or polyacetylene. It has been proved 
that for a closed ring of N atoms, with Hamiltonian 
as in ID), the minimum energy configuration is precisely 
the one in which there are two alternating bond lengths 
h and h, provided N = 2 mod 4 (TV = 6, 10, 14, . . .). It is 
possible that h — h, as is the case in benzene, but noth- 
ing more complicated than dimerization can ever occur 
in this model, for any F-function. When N = mod 
4 more complicated things can occur when F is not a 
quadratic, e.g., for A'^ = 4 the minimum configuration 
can be a trapezoid [l]. Nevertheless, as N goes to infin- 
ity the dimerized energy is always asymptotically exact. 
The general case is analyzed in [G] in detail and the stated 
conclusion is even shown to apply to the so-called 'Spin- 
Peierls model'. 

The method of proof for ID does not generally extend 
to higher dimensions. It does not extend to the square 
lattice unless a physically unreachable magnetic field is 
imposed but, as we discover here, it does extend to 
the hexagonal lattice with zero magnetic field. The rea- 
son for this difference, ultimately, is that the energy min- 
imizing magnetic fiux for the square lattice is non-zero, 
whereas for the hexagonal lattice it is zero, as assumed 
in H(T) above 0. 

Naturally, we cannot expect perfect periodicity in a 
finite lattice, because there will always be edge effects. 
This remark motivates the following theorem, whose 
proof we will outline here. 

Theorem. — For any triplet of hopping matrix el- 
ements (t(i),t(2)^^(3)) e(t(i),i(2)^i(3)) denote the 

ground state energy per atom in the thermodynamic limit 
for the Hubbard model ^ with these t 's arranged as in 
Fig. 1(a). Let e be the minumum of e{t^^\t^'^\t^^^) with 
respect to these t 's. Then the ground state energy E(T) 
of ([T]) with arbitrary t's on a large, finite lattice A, with 
|A| sites, satisfies 

E{T)>\A\e + o{\K\), (2) 

where o(|A|) < const |A|"'^/^ In |A|. The same conclusion 
holds for the free energy at finite temperature. 

Strictly speaking, we prove that there is a minimizer 
(to leading order) with the structure of Fig. 1(a), but 
we do not rigorously prove uniqueness of this structure. 
Nevertheless the uniqueness is evident for two reasons: 
One is that our proof proceeds by showing that repeated 
reflection of a single line cannot raise the energy, but in 
reality it surely lowers it. The second is that the struc- 
ture made this way out of repeated reflections still has 
the hopping matrix elements that were in the given row 
before it was reflected. At this point we could lower the 
energy still further by reevaluating the optimum i- values. 
These will surely change, and the energy will be lower 
because the equation for the optimum values, namely. 



F'{txy) — 2(c|.Cy) = dE/dtxy = will likely no longer be 
satisfied with the old t values. For example, the ALT 
structure in Fig. 1(b) morphs into the Fig. 1(a) structure 
of the Kekule type. If the t value on the heavy bond did 
not change, then the expectation value of cJ,Cy along the 
distinguished bond would have to be the same as for the 
undistinguished bond - which it is not. 

Wc recall two facts about the Hubbard model H]), 
which hold even for arbitrary T on any finite, bipartite 
lattice A with an even number of sites. 

(1) Among the absolute ground states there is one with 
N — |A|=number of sites. If U ^ 0, this state is unique 
6] . This means that we do not have to worry about con- 
straining the particle number to be equal to the number 
of lattice sites, |A|. This constraint is automatic. 

(2) Am ong the ground states there is one with total 
spin zero |10|. 

The main tool we are going to use is reflection positiv- 
ity which originated in quantum field theory [ll[ . It was 
later used in statistical mechanics (l^ - [l6l | to prove long 
range order of systems with continuous symmetry. It was 
also important for static situations like the solution of the 
flux phase problem 0]. It is in the latter sense that we 
apply it here. Our theorem above can be rephrased by 
saying that the lowest energy occurs for a lattice config- 
uration that preserves the three-fold reflection symmetry 
of the hexagonal lattice, as shown in Fig. 1(a). 

Let A be piece of a hexagonal lattice that is reflec- 
tion symmetric about a line I that passes perpendicu- 
larly through bonds. The two halves of the lattice will 
be called A^ and A^. Let Tl and Tr denote the hopping 
elements on the left and on the right side, respectively, 
and let Tm be those corresponding to bonds cut by /. 
Let T£ be the reflection of Tl to the right and, simi- 
larly, let be the reflection of Tr to the left. Then, 
symbolically, T = {Tl^Tm^Tr) and we want to consider 
= {Tl,Tm,TI) and T^ = (T^.Tm^Tr). The corre- 
spending ground state energies are related by the basic 
inequality, whose proof is in [9, ^ [l3| ; 

E(T) > \E(T^) + \E{T^) . (3) 

Our way forward is now obvious. Think of an in- 
finite hexagonal lattice. There are three directions in 
space for which we can draw lines in these directions 
that cut lattice bonds perpendicularly (but do not pass 
through sites). For each of these three directions there 
are infinitely many parallel lines with these property. We 
would like to use reflection positivity for each of them to 
say that if the hopping matrix elements are not as shown 
in Fig. 1(a), then we can lower the energy by reflecting 
one half of the lattice onto the other half and keeping the 
better of the two choices according to inequality ([3]) . In 
short, if the optimum configuration did not have the sym- 
metry of Fig. 1(a) the energy could be lowered, thereby 
giving a contradiction. 



The physics is now clear. The mathematical problem 
is to make rigorous sense of the lattice reflections by con- 
sidering finite lattices and passing to the thermodynamic 
limit in an appropriate way so that at each stage the er- 
ror we make is of lower order than the area |A|. This is 
accomplished by imposing periodic boundary conditions 
on the finite lattice, but it is much more complicated, 
geometrically, than a similar problem for the square lat- 
tice (with an unrealistic magnetic field flux). There, one 
can impose periodic boundary conditions in two orthog- 
onal directions simultaneously. This is impossible for the 
hexagonal lattice since the three reflection directions are 
not orthogonal to each other, and we can impose periodic 
boundary conditions in only one direction at a time; see 
[l6| for a discussion of this. 

Nevertheless the problem can be solved and we sketch 
the solution here for the ground state energy. The exten- 
sion to finite temperature follows by similar arguments 
using the reflection positivity established in [9|, Il3l4l5l | . 

Proof sketch of the Theorem. — Step 0. The con- 
clusion shown in Fig. 1(a) can be deduced from a state- 
ment about individual bonds. Any bond, &, in our lattice 
(except for those at the boundary) has four other bonds 
connected to it, as shown in Fig. 2(b). Two of these, to- 
gether with b, form part of a hexagon and the other two, 
together with 6, form part of a neighboring hexagon. We 
say that b is happy if the first two bonds have equal t- 
values and also the second two have equal t-values (pos- 
sibly different values) . The i- value of b itself is irrelevant 
for this definition. (A similar definition applies to bonds 
at the edge, which have fewer than four bonds attached 
to them.) 

Our goal is to show that every bond in a minimizing 
configuration is happy, for this will imply the theorem. 
To be precise we shall show that all except lA|^/^ln|A 
of the bonds are happy, but this is good enough for the 
thermodynamic limit. 

Step 1. We consider an increasing sequence of A's each 
of which is a rectangular piece of the hexagonal lattice, 
with even side lengths of order L. Approximately one 
third of the bonds are oriented in the y-direction. We 
periodize A in the y-direction by connecting the sites on 
the top row to the corresponding sites in the bottom row 
to form a new row of hexagons, and call the new lattice 
A. This leads to a cylinder whose axis is parallel to the 
X-axis. The new hopping matrix T consists of the old 
hopping matrix T and some new positive elements txy 
connecting the top to the bottom. The new Hubbard 
Hamiltonian H{T) is defined as in ([T]) except that the 
(positive) ixy's which connect top to bottom are inserted 
with a minus sign. With this convention, every hexagon 
on the cylinder has flux zero and any circuit around the 
cylinder that encircles the cylinder axis has length mod 
4 and has flux tt (The flux of a circuit is the argu- 
ment of the product of the Vs along this circuit.) The 
importance of these fluxes lies in the following fact: Re- 






FIG. 2: Some definitions: (a) The top left figure is a single 
row and the bottom left figure is a reflected 6 row. They 
are joined by overlapping the common vertical bonds, (b) A 
happy bond (heavy line) has the values of the four adjacent 
bonds equal in pairs, as shown. 



flection positivity requires dividing the lattice into two 
pieces and this requires cutting bonds, which means cut- 
ting closed circuits. In order to have reflection positivity 
every closed circuit that is cut must have flux tt if the 
length of the circuit is mod 4 and flux if it is 2 mod 4. 
This is shown in [13] . The length of a hexagon circuit 
is 6 and the length of a circuit going around the cylinder 
is always mod 4. 

The insertion of additional bonds changes the energy 
from E{T) to £^(7'), the ground state energy of H{f ), by 
at most 0{L), i.e., order 1/L per atom. 

We think of this cylinder as composed of rows, parallel 
to the axis, as shown in Fig. 2(a), together with their 
reflections, which we shall call 6'-rows. Thus A is a se- 
quence of rows interspersed with 9 rows. Each row and 9 
row has its prescribed hopping matrix elements. Let us 
denote the rows as we go around the cylinder by Ai, 6A2, 
A3, 9A4, Am-1, 9Am (with M of order L). From 
these we can create M new configurations Ti, . . ., Tm 
of matrix elements of which the j-th one has the rows 
Aj , 9 A-j , A-j , 9A.j, .... It is a well known consequence of 
reflection positivity ^ (called the chessboard estimate 
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(4) 



from which we see that one of the Tj configurations is 
energetically at least as good as the original T configura- 
tion. In a Tj configuration every vertical bond is happy. 

At this point we have established a T configuration 
for which all the vertical bonds are happy. We have no 
knowledge of the happiness of the zig-zag bonds between 
the vertical bonds. Our next goal is to make 50% of the 
remaining zig-zag bonds happy without destroying the 
happiness of very many of the vertical bonds. This is the 
difficult step. 

Step 2. We resist the temptation to cut the cylinder 
along the line by which it was originally formed. Instead, 
we proceed as follows, and we urge the reader to follow 
this construction by forming a cylinder from a piece of 
paper and tape. We start at one edge of the cylinder 
and, with a pair of scissors we cut the paper along a 
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FIG. 3: (a) The heavy hnes in the left figure show the more 
than 50% happy bonds after reflecting in orientations 1 and 
2. (b) shows the more than 75% happy bonds after reflecting 
in orientations 1, 2 and again 1. 

line that is 30° from the axis (so 60° froni the vertical 
direction) . The line should not go through the atoms and 
it should cut an even number of bonds. In this way we 
cut a spiral that eventually reaches the other edge of the 
cylinder and we lay the paper out flat on the table. It 
will be noted that there are two ways to tape this lozenge 
piece of paper together to form a cylinder. One way is 
just to put it back together in the way we cut it. The 
second way, which is the one we will use, is to connect 
the lozenge to itself with a shift in the x-direction. We 
pay a surface energy to accomplish this. The old vertical 
bonds that were parallel to the y-axis now lie at an angle 
of 60° to the axis of the new cylinder (which we again 
assume to be parallel to the x-axis), and a different set 
of bonds now lies parallel to the y-axis. 

Again we focus our attention on a single row of the 
new cylinder, as shown in Fig. 2(a). The happiness es- 
tablished in Step 1 means that every second zig-zag bond 
of the new row is happy. We now repeat the reflection 
positivity argument. Inequality ^ allows us to repeat 
one row periodically. This repetition preserves the hap- 
piness of every second bond on the zig-zag line and makes 
every vertical (in the new orientation) bond happy. The 
resulting network of happy bonds consists of chains as 
shown by heavy lines in Fig. 3(a). 

Step 3. We repeat Step 2 by cutting a spiral and re- 
tying the cylinder as it was initially in Step 1. Since we 
have cut the chains from Step 2 at an angle of 30° relative 
to the cylinder axis, we observe that every row parallel to 
the cylinder axis has a zig-zag line with three out of every 
four bonds happy. Reflection positivity dU allows us to 
repeat one such row endlessly. All the resulting vertical 
bonds are happy and out of every sequence of four zig- 
zag bonds three are happy. The configuration of happy 
bonds again consists of chains but now these chains are 
thicker; on every row the chain has three zig-zag bonds 
instead of one. This is shown in Fig. 3(b). 

Step 4- We cut and paste again and use reflection 
positivity to achieve happiness for seven out of eight zig- 
zag bonds. The chains are now seven bonds thick. 



After k steps all the vertical bonds and 2*^ — 1 out of 
every consecutive 2*^ zig-zag bonds are happy. We have 
paid a surface energy k times and we also might have 
incorrect bonds along the sequence of cut lines, but each 
of these corrections can change the energy at most by a 
quantity of order L. If we take k hi2 L, we find that all 
except possibly order kL bonds are happy, and we have 
paid an energy error of at most 0{L ln2 L). This is what 
we wanted to show. 

Conclusion. The Hubbard model for graphene on the 
hexagonal lattice has a special property called reflection 
positivity that the square lattice does not have. With the 
aid of this property we are able to limit the kind of lattice 
distortions allowed by this model. Only distortions that 
preserve the three-fold reflection symmetry of the lattice 
are possible in the thermodynamic limit. 
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